Rayleigh- Taylor instability of crystallization waves at the superfluid-solid ^He interface 
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At the superfluid-soUd *He interface there exist crystallization waves having much in common 
with gravitational-capillary waves at the interface between two normal fluids. The Rayleigh- Taylor 
instability is an instability of the interface which can be realized when the lighter fluid is pro- 
pelling the heavier one. We investigate here the analogues of the Rayleigh- Taylor instability for the 
superfluid-solid ''He interface. In the case of a uniformly accelerated interface the instability occurs 
only for a growing solid phase when the magnitude of the acceleration exceeds some critical value 
independent of the surface stiffness. For the Richtmyer-Meshkov limiting case of an impulsively ac- 
celerated interface, the onset of instability does not depend on the sign of the interface acceleration. 
In both cases the effect of crystallization wave damping is to reduce the perturbation growth-rate 
of the Taylor unstable interface. 

PACS numbers: 68.03.Kn, 67.80.-s, 52.35.Py 



INTRODUCTION 

The Rayleigh- Taylor instability is a fingering instabil- 
ity of an interface between two fluids of different den- 
sities. It takes place when the heavier fluid is deceler- 
ated by the lighter fluid or, in other words, density and 
pressure gradients have opposite directions [l|, The 
Rayleigh- Taylor instability 3] occurs in numerous phys- 
ical and technological situations, e.g., gravity-driven in- 
stability of a heavier fluid atop a lighter one and inertial 
confinement fusion. In essence, the Rayleigh- Taylor in- 
stability is the first step in a fluid-mixing mechanism, the 
step eventually leading via formation of bubbles, spikes, 
and curtains to the turbulent regime of fluid mixing. 
Along with the Kelvin-Helmholtz criterion for tangen- 
tial discontinuities at the interface between two normal 
fluids, the Rayleigh- Taylor criterion is among the most 
generic principles in the complicated subject of interface 
instability. 

For the interfaces of superfluid ''He, some of the hydro- 
dynamic instabilities have been found as well. First, we 
mention the Faraday instability which denotes the phe- 
nomenon of the parametric excitation of standing waves 
on the free surface of a fluid. The flat shape of the sur- 
face becomes unstable with a periodic modulation of the 
acceleration of gravity. Recently 0] , generation of Fara- 
day standing waves on the free surface of "^He has been 
realized in the experimental cell subjected to sinusoidal 
vibration in the vertical direction. According to 0, @ , it 
is also possible to generate a dense fog of helium droplets 
by driving the capillary waves on a superfluid ^He surface 
unstable with an intense ultrasonic beam from a piezo- 
electric transducer under the surface. There have been 
observed electrohydrodynamic interface instabilities due 
to charges trapped at the surfaces and interfaces of vari- 
ous condensed helium phases [3l . The shear flow between 
the superfluid A and B phases of '^He can result in the 
Kelvin-Helmholtz interface instability Q. 

Dynamics of the superfluid-solid ''He interface due to 



sufficiently fast processes of crystallization and melting 
resembles much that of the free surface of a fluid. In 
particular, as was predicted by Andreev and Parshin, 
the crystal in contact with its liquid phase can support 
wave-like processes of crystallization and melting, see re- 
view [9J . From the dynamical point of view such weakly 
damping crystallization waves are an immediate coun- 
terpart of the well-known gravitational-capillary waves 
at the vapor-liquid interfaces. 

A series of mechanical and hydrodynamical instabili- 
ties has been predicted and observed for the superfluid- 
solid ''He interface. We mention the Grinfeld instability 
under uniaxial stress of a solid. Warping of the flat in- 
terface occurs at some threshold stress when the release 
of elastic energy exceeds the loss of the surface energy 
Like normal fluids, the steady flow of a superfluid 
in the direction tangential to the interface can result in 
the Kelvin-Helmholtz instability. As the flow exceeds a 
threshold magnitude, crystallization waves appear at the 
superfluid-solid "^He interface [HI, ■ The phenomenon 
has qualitatively been observed as a distortion of the 
crystal surface in the fluid jet [l^]. An analogy with 
generating sea waves by wind is fully appropriate here. 

To date, the Rayleigh- Taylor phenomena have exten- 
sively been studied in normal fluids, but not much study 
has been made in superfluids or quantum solids. The 
classical Rayleigh- Taylor instability of the superfluid- 
solid ''He interface in the field of gravity is observed by 
Demaria, Lewellen, and Dahm IJ) . In these experiments 
a cell in which the solid and liquid phases occupy initially 
the lower and upper halves, respectively, is inverted me- 
chanically by 180°. After inversion a single finger of the 
liquid phase ascends at the centre of a cell, and the solid 
phase descends along the walls. On the other hand, in 
experiment [l5| the flat shape of the interface remained 
stable for a ''He crystal grown at the needle point with its 
lower facet under favorable condition for developing the 
gravity-driven instability. A difference in the observa- 
tions can be associated with the following reasons. The 
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FIG. 1: The growth of a ''He crystal at 0.47 K and initial 
overpressure 5.2 mbar. The left frame corresponds to 0.19 ms 
after nucleation. The right frame is taken 80 ms later when 
the net pressure is already close to the melting pressure. The 
vertical size of the frames is 2.4 mm. 



requirement for the interface instability as well as the 
initial stage of fingering process, as is shown in [l6| . is 
sensitive to the state of a crystal facet, rough or smooth, 
and to the size of a facet. 



Recently [17| , it has been demonstrated that the Fara- 
day instability is also inherent in the superfluid-solid ^He 
interface. Crystallization waves at the horizontal inter- 
face between superfluid and solid *He are generated by 
a periodic vibration of an experimental cell in the ver- 
tical direction. In accordance with expectation 18| the 



amplitude of the waves excited at one-half of the driv- 
ing frequency decreases for higher temperatures due to 
reduction of the interface growth coefficient describing 
dissipative properties of the interface. From general ar- 
guments the Faraday instability can be viewed as a par- 
ticular case of the Taylor instability for the periodically 
driven interface. 

On the other side, the spectrum of crystallization 
waves remains invariable for the steady flow of a super- 
fluid in the direction normal to the interface 12| . At first 
sight this implies that the growth of a solid should not 
bring the superfluid-solid interface to instability. How- 
ever, in the experiments on free growth of a ^He crystal 
initiated at the need le p oint immersed into the overpres- 
surized liquid bulk [l^ EO] one can observe a destruc- 
tion of the regular shape of the crystal triggered under 
sufficiently large overpressure exceeding about 6 mbar. 
Immediately after nucleation the crystal seed has a clear 
hexagonal prism-like shape with slight ripples. Soon af- 
terwards by 0.1-0.2 ms the shape of the crystal becomes 
round with a highly irregular outline. Far later by 100- 
150 ms, as the net overpressure vanishes and the pressure 
in the cell becomes phase-equilibrium, the shape of the 
grown crystal relaxes to a regular hexagonal prism. We 
put two images of a crystal in Fig. [TJ 

For overpressures higher than 6 mbar, see Fig. [21 we 
discover more exotic patterns by the same time interval 
0.1-0.2 ms after the voltage pulse which triggers nucle- 
ation. The interfacial irregularities become more pro- 
nounced and acquire a mushroom-like shape. Appar- 
ently, the fluid moves into the crystal between the neigh- 
boring spikes, resulting eventually in generation of liquid 
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FIG. 2: The deviation of the pressure from equilibrium vs 
time during crystal growth at 0.48 K and initial overpressure 
7 mbar. The left insert shows a crystal at 0.64 ms after nu- 
cleation. The time at the right insert is 16 ms. The vertical 
size of the frames is 2.4 mm. 



bubbles inside the crystal. 

In addition, in Fig. [2] we show the behavior of the pres- 
sure in the course of the crystal growth. After nucleation 
of a crystal seed the overpressure in the cell drops and 
then gradually vanishes, oscillating around zero value 
corresponding to the equilibrium pressure. The pres- 
sure oscillations are due, in the first turn, to the finite- 
ness of the experimental volume and the finiteness of 
sound velocity. The point is that the appearance of a 
solid seed and its next growth are accompanied by varia- 
tions of the density and the volume in which the density 
changes. The change of the seed volume results in ex- 
citing and emitting the sound waves which propagate in 
the direction to the container walls with the next back- 
ward reflection from the walls to the solid seed. E.g., 
the first- and second-sound emissions with an expanding 
■^He-concentrated drop in a superfluid ■^He-'*He mixture 
has been analyzed in Ref. 2]J. The excitation and emis- 
sion of sound become more and more effective as the 
interface rate and acceleration increase. Eventually, we 
obtain acoustic damping oscillations of the liquid inside 
the cell [2^ . In some sense the damping of the pressure 
oscillations represents a quality factor of a liquid / solid or 
melting/freezing resonator. In the lack of any dissipative 
processes in the system the pressure oscillations will last 
infinitely long. 

The pressure in the cell becomes equilibrium pressure 
in 2 ms and the driving force vanishes. Finally at 16 ms, 
the crystal relaxed and acquired the regular hexagonal 
shape. However, we still observe the liquid bubbles em- 
bedded into the crystal bulk. 

If we roughly estimate the velocity which the interface 
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should acquire by the time of the overpressure release and 
formation of the irregular outline, we will find a rather 
high magnitude of several meters per second. The cor- 
responding acceleration which provides such increment 
of velocity should be about lO^g, g being acceleration of 
gravity. On the whole, large acceleration for a short time 
shows evidence in favor of a shock-driven character of the 
crystal growth. It is interesting to note that the irregu- 
lar patterns observed are similar in appearance to those 
obtained in studies [2§| of a shock-accelerated boundary 
between two gases of difi'erent densities. The typical time 
of the pattern formation was of the same order of several 
tenths of millisecond. 

In the present work we develop a linear theory on the 
Taylor instability of an arbitrarily accelerated boundary 
between the superfluid and solid ^He phases. In essence, 
we derive an equation which the small interfacial pertur- 
bations obey. We consider three typical cases of the in- 
terface acceleration: constant, shock, and periodic. The 
plane and spherical interface geometries are analyzed. 



PLANE GEOMETRY AND THE INTERFACE 
GROWTH KINETICS 

Let us assume the interface is parallel to the x-y plane, 
with vertical position z = L{t) which moves at the rate 
V = L{t). The upper half-space z > L{t) is occupied 
with the liquid phase and the solid phase occupies the 
lower one z < L{t). Below we will consider the stability 
of the moving interface with respect to its small perturba- 
tions C,{x,y, t) from the flat shape. Thus, the coordinate 
Z[x, y, t) = L{t) + C,{x, y, t) gives the vertical position of 
the perturbed interface evolving in time. 

To discover the effect of nonuniform motion on the sta- 
bility of the interfacial shape, we make a number of sim- 
plifying assumptions which do not affect the main point 
of the phenomenon. The validity and criteria of appli- 
cability for the assumptions to be made below can be 
found in the papers on the kinetic interface coefficients 
and crystallization waves (23 - [27| . 

So, in the superfluid we employ the usual two-fluid 
equations without dissipation. In addition, we assume 
that the growth rate of a solid V{t) is always small com- 
pared with the velocity of the first or second sound. This 
is an ordinary experimental situation. Thus we treat the 
hydrodynamics of the superfluid in the approximation of 
incompressible liquid and the constancy of the entropy 
density. In this case [ll] the velocities and of the 
normal and superfluid motions can be described in terms 
of gradient of velocity potentials 0„ and (ps which satisfy 
V^0„ — and V^t/i^ = 0, respectively. Since flnally we 
will discuss only the linearized equations in the perturba- 
tion it is convenient to consider a single Fourier mode 
of the perturbation ^ = C,q{t) exp(iqr) with wave vector 
Q = {lx, Qy) parallel to the boundary. The solutions of 



can be represented as 

4is = Us{t)z + As{t) exp{iqr - qz), 
(t>n = Un{t)z + An{t) exp(zqr - qz), 



(1) 



where r = {x,y) and the velocities Us{t) and Un{t) de- 
scribe the undisturbed motion in the superfluid. The per- 
turbation amplitudes As{t) and An(t) are assumed to be 
linear in Q (t) and will be determined later from the corre- 
sponding boundary conditions at the interface. The pres- 
sure in the superfluid is a sum of pressures P = Pn + Ps 
and 



Ps =Psoo- Ps{k + (V(/).)V2) - Psgz, 
Pn = Pnoo~ Pn{(l>n + [V(j)nf /l) ~ PnQZ. 



(2) 



The index "c»" refers to the values taken away at infinity. 

Unlike previous considerations which are also 

linear in the interfacial perturbations, we have to retain 
the quadratic terms in the superfluid and normal veloc- 
ities Vs.n = V0s,n on account of nonzero value Us.n(t) 
and product like UsAs or UnAn- 

The mass continuity across the boundary gives at z = 
Z{x,y, t) 



ju = PnVn V + PsVs i/ = (P - p')^ ) 



(3) 



where Z is the velocity of the boundary, i> is the unit 
vector normal to the boundary, and ji, is the mass current 
normal to the boundary. The densities p„, Ps are the 
normal and superfluid densities, p = p„ -I- Ps, and p' is 
the density of the solid phase. The normal components of 
velocities can be approximated by « ti„ 2 and i/ ~ 
Vs z ■ We also believe that there is no motion in the solid 
phase, i.e., w' — 0. 

To further simplifications, we suppose the normal com- 
ponent sticks to the interface like a viscid fluid. Also, this 
implies the Kapitza resistance to be infinite [25| . In ad- 
dition, we disregard any excitations, e.g. phonons, in the 
solid. So, we put at the boundary 

Vn.^Z. (4) 

Using Eqs. ([T]), and ([4]), we can determine the un- 
known amplitudes in ([T]) and then the velocity fields Vs 
and Vn in the liquid. For the unperturbed motion, it is 
obvious 



P - Ps 



V and Un — V . 



(5) 



The amplitudes An (t) and As (t) are given approximately 
by 



Asit) 



_ P' - Ps Cqjt) „L{t) 



^„(t) = _M)e«^(*) 
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As a result, we can also calculate the pressure field ([2]) in 
the liquid. The boundary conditions ([3]) and (HJ in com- 
bination with the obvious relations (O determine unam- 
biguously the inertial properties of the interface described 
in terms of the effective density [25] 

Pel = Pn + {p' - PsY IPs 
in the sense that 

To proceed further, we adopt the most simplifying as- 
sumptions [255 to describe the solid and its boundary. 
The solid is assumed to be always unstressed and all 
possible shearing components Ui^k of the stress tensor 
Gik are neglected. In other words, the stress tensor is 
isotropic, i.e., aik = —P'Sik, and we can define 'pressure' 
according to P' — —(Tu/S [27]. Then, from the formal 
point of view, the solid can be described as a liquid un- 
der pressure equal to P' . 

The next boundary condition stems from the conti- 
nuity of the momentum flux density across the inter- 
face. The momentum flux density in the superfluid 
[2^ reads PSik + PnV„iVnk + PsVsiVsk- Then, we take 
'^ik^k = —P'vi into account, assume the small curvature 
of the interface z = Z{x, y, t) = L(t) + C,(x^ y, t), and use 
a frame that refers to the boundary 

P + Pn{Vn - Zvf + Ps{Vs - Zvf - [P' + p' Z^) = 

= -l.kdZ-'/dndrk = l^kd^'/dndrk . (6) 

Neglecting the quadratic terms in velocities gives the 
usual Laplace condition of mechanical equilibrium across 
the interface |25|. Here 7zfc(^, y) = a5ik+ da^ / dipidipk is 
the surface stiffness tensor [§, llO, [l^ expressed in terms 
of surface tension a = a{9, (p) depending on the angles 
between the crystalline orientation and the normal to the 
surface. 

Let us turn now to the last boundary condition. It is a 
reasonable assumption that any motion of the interface 
accompanied also by the melting and growth of a solid 
will dissipate a certain amount of energy. Thus a finite 
velocity of the interface should produce some imbalance 
in the chemical potential difference p — p' between the 
liquid and solid. The routine in various theories of the 
interfacial dynamics is an introduction of the so-called 
growth coefficient K which relates the interface growth 
rate with the difference in chemical potentials across the 
interface 24, |25[. Because of Un{t) 7^ and Us(t) ^ 
we again have to take into account the squares of veloci- 
ties which are always omitted in the linear perturbation 
theory of the interface being initially at rest. So, at the 
boundary we employ an effective relation 



K 



{Vs - Zv) 




(7) 



where p and p! are the chemical potentials of the liq- 
uid and solid per unit mass. In a wide sense the growth 
coefficient here is a certain combination of all Onsager co- 
efficients and the kinetic coefficients describing the near- 
surface dissipative processes. In general, the growth co- 
efficient K can depend on the temperature as well as on 
the wave vector q. Usually, in the ballistic regime, when 
the mean free path I of excitations is large, the growth 
coefficient is independent of wave vector. In the opposite 
hydrodynamic limit ^ 1 the growth coefficient may 
defend on the wave vector approximately as XjK ^ ql 

Lastly, we need an expression for the chemical poten- 
tial difference. As usual, the reference point is the melt- 
ing pressure Pc at which the chemical potentials p and 
p' coincide and the liquid-solid transition takes place. 
We take the necessary formulae for the superfluid from 
Ref. [2^. After expanding chemical potentials in the 
vicinity of the melting pressure, we obtain 



= a{T - Too) 



p-p = 

P-Pc Pn {Vn - Vs)^ P' - Pc 

p p 2 p' 



crp 



Pn 



Ps 



where a is the entropy and the quantities with index "cx)" 
stand for the magnitudes taken far from the interface. 

Now we are in position to find the equations which 
the interface dynamics obeys. Knowing velocity poten- 
tials (j)„ and 05 expressed via and perturbation 
Cq{t)j we can calculate the normal and superfluid veloc- 
ities, pressure, and chemical potential difference. Next, 
we insert the quantities calculated at the interface into 
the boundary conditions ([6]) and ([7]) and eliminate the 
pressure P' . As a result of some algebraic formula ma- 
nipulation linear in (^q, we obtain an equation consisting 
of the -independent component and the other one linear 
in (^q. The first component gives an equation 



K 



[AP - pgL)] 



Pel 



y2 

VL+- 



which describes the undisturbed motion of the flat in- 
terface and relates overpressure AP{t) = Poo{t) — Pc to 
V{t) = L(t) in a complicated manner in order to support 
the necessary behavior of the growth rate. This equation 
does not have much interest for us. 

The other equation obtained is the most significant 
one. It represents the equation for the linear dynamics 
of the interface perturbation C = C,q{t) exp(jqr) when the 
interface is subjected to an arbitrary driving acceleration 
V{t) 



Pel 



- + ^ 4 + hkmk + (p' - P)9 - PeiV{t)\ Cg = 0. (8) 
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The new aspect of the equation derived is an additional 
term with the interface acceleration V{t). As is expected, 
the uniform motion of the interface at V{t) — const does 
not influence the character of small interfacial oscilla- 
tions. For V{t) — 0, Eq. amounts to the known re- 
lation determining the spectrum of crystallization waves 
when the melting-crystallization processes are balanced 
and the interface position in average is invariable, i.e., 
L{t) = const [i,[2Sl2i- 

Undoubtedly, more realistic and complicated models of 
the superfluid-crystal "^He interface will improve the mag- 
nitudes of the effective interface density, effective growth 
coefficient, and surface stiffness. However, we believe 
that the structure of Eq. ([8]) is generic and holds. 



INTERFACE INSTABILITIES 

Equation (|S]) can have unstable solutions depending 
strongly on the acceleration history of the interface. 
First, let us consider the stability of the plane interface 
with respect to small perturbations Cg ~ exp{X{q) t) for 
the uniformly accelerated growth of a crystal. The root 
with ReA(g) > means the interface instability, i.e., ini- 
tially small-amplitude perturbations of wavelength 27r/g 
will grow exponentially in time. For the acceleration ex- 
ceeding the threshold Vc — g{p' — p)/ PcU the interfacial 
perturbation will increase for the wave vectors satisfying 
79c < PcfV^ - [p' - p)g- For brevity, we put jikqiQk = 79^ 
where 7 = jiknin^ and n — q/q is unit vector in the 
direction of perturbation. Note that the threshold ac- 
celeration Vc does not depend on the growth coefficient, 
i.e., on the dissipative properties of the interface, and is 
positive. The latter corresponds to the case when the 
interfacial acceleration V is directed to the superfluid. 
Thus, the Taylor instability due to non-uniform growth 
of the plane interface appears only for the accelerated 
growth of a solid. 

The surface stiffness stabilizes the region of long wave 
perturbations and establishes the most unstable wave- 
length having the fastest exponential growth. The value 
qo corresponding to the maximum magnitude ReA(g) 
gives the shortest time scale for the development of the 
instability which will be characterized by the wavelength 
27r/(7o- The value go can be found from the equation 



approximately equal to 



% - 3 



V2 p' 
3 



3/2 
go 

K 



where q^. is the value related to the upper bounds of in- 
stability according to 79^ = p^iV — g{p' — p). For large 
magnitudes of the growth coefficient or large acceleration 
in the case K'^V ^ p'"^ l^p\^ , the values qo and X(qo) are 



- 

\3V3 Pcf 



i6VS-fPeiK^qcV/^ 

.3 \ 1/2 



Pe{K 



The values go and Aq depend on 7, thus implying 
anisotropic and complicated possible surface patterns. 
In the opposite limit K*V <C p''^ llP% one has roughly 



go = 

Ao = — 7'7c 
P 



IPol 

K 



1/3 



-(pef^-(p'-p).g) 



It is interesting that, though the spatial scale g,^^ is sen- 
sitive to the surface stiffness and its anisotropy, the time 
of developing the instability becomes independent of the 
surface stiffness and its anisotropy. 

On the whole, the values go and Ao decrease as the 
kinetic growth coefficient reduces or dissipation with the 
interface enhances. From the experimental point of view 
this may require a crystal surface of sufficiently large sizes 
d > 27r/go and a large time to support the accelerated 
growth t > 1/ Aq in order to realize an interfacial insta- 
bility with a uniformly accelerated growth. 

In Ref. Is*] it has been shown that periodic modulation 
of the gravitational constant as g{t) — g{l + 2ecos2ujt) 
can result in parametric excitation of crystallization 
waves at the stationary flat interface corresponding to 
V(t) = in our case. Admitting some analogy between 
gravity and noninertial frame, one can expect a possibil- 
ity of exciting crystallization waves with a periodic driv- 
ing, e.g., V{t) — Gcos2ujt, even with the lack of gravity. 
To demonstrate this, let us rewrite Eq. ([5]) in the form 



Q + rgQ+LoUq){l~V{t)/g)Q = 0, 



(9) 



where we have introduced the damping coefficient Tq, 
the frequency a;o(g) of crystallization waves in the lack 
of damping, and the scaled accelerating amplitude g 

^liq) = {iikqiqk + {p' - p)g)q/pcu = — , 

Pcf 



Pcf 



likqiqk 
ip' - p)9 



Then, if a periodic process of melting and crystallization 
is realized in experiment so that the interface could oscil- 
late around some average position at frequency 2u and 
amplitude G/(2w)^, equation ^ with label 2e = —G/g 
transforms to a Mathieu-type equation 

'(q + TqC,q+ujl{q){l + 2e COS 2a;i) Cg ^ 0, 

which is identical to that analyzed in [l8|. The predic- 
tions which follow are well known and we refer the readers 
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to papers [17|, [18| for details. Note only that in the free 
crystal growth experiments 15, [23| the pressure in the 



have 



cell drops drastically down after nucleating a solid seed 
with the subsequent transition to the damped oscillations 
around the melting pressure. 

Crystallization waves at the superfluid-solid ^He inter- 
face can also be generated with the Richtmyer-Meshkov 
mechanism when the interface is subjected to an impul- 
sive acceleration, i.e., V{t) ~ V{0)S{t). In ordinary fluids 
and gases, for this purpose, the passage of a shock wave 
across the interface is commonly used. Unlike the Taylor 
case of constant acceleration when the perturbation am- 
plitude in the linear regime grows exponentially in time, 
the initial stage of interface instability in the Richtmyer- 
Meshkov case of shock acceleration |29[ is characterized 
by a linear growth of the perturbation amplitude in time. 
The Richtmyer-Meshkov instability is independent of the 
direction of acceleration in contrast to the Taylor one. 
The late time stages of both instabilities may demon- 
strate a formal resemblance, showing bubble and spike 
morphology. 

According to Eq. ([8]), the growth rate of the Richtmyer- 
Meskov unstable interface in the linear regime can ap- 
proximately be described in the overdamped regime as 



sinh('i^r2/4-c.2(q) 



(10) 



and in the weakly damped regime as 
sinft^u;^,{q)~rl/A 



^ul{q)~Tl/A 



C,(a;,0), uJo{q)>Tj2. 



Here Cq(x,0) — Cq(0) cosga; is the initial perturbation 
amplitude at the interface immediately after the impul- 
sive acceleration and V^(0) is an increment in the interface 
velocity due to acceleration of the boundary. At the lin- 
ear stage the amplitude of crests and troughs is approx- 
imately the same. The shape of the crests and troughs 
is similar and the interface remains approximately sinu- 
soidal. Staying in the linear regime, we see from ()10|) that 
the finite damping of crystallization waves is a stabiliz- 
ing factor and can limit the growth of the perturbation 
amplitude. 

Next, we consider two extreme cases of wol?) ~ and 
— 0. For a;o(g) = which corresponds mainly to the 
fiat interface perturbations with small wave vectors g, we 



= yK ViO) i^l - e-^ j (.{x , 0), ^.0(9) = . 

Here it is worthwhile to note that the total amplification 
factor is proportional to the interface growth coefficient 
K and proves to be independent of the wavelength 2n/q 
of an initially sinusoidal perturbation. On the contrary, 
the time for ceasing the growth of perturbation amplitude 
is scale-dependent. The larger the wavelength, the longer 
the decay time. A special feature of crystallization waves 
in ^He is that the growth coefficient K, and thus Fg, is 
strongly temperature-dependent. 

For the opposite case when the crystallization wave 
damping is absent Tq = 0, we arrive at the persistent 
oscillations of the initial interface perturbation 



i{tuJo{q)) 
wo(f?) 



Cg{x,0), F, =0. 



The maximum amplitude of oscillations is governed 
by the crystallization frequency wo(<7). Since roughly 
wo(?) ^ <Z'^^^, the large curvature of the initial interface 
perturbation can also be a factor that limits the infinite 
growth of perturbation amplitudes. 

In ^He the Richtmyer-Meshkov instability can be stud- 
ied with a sound wave pulse hitting the interface in the 
normal direction. Recently [30], it has experimentally 
been demonstrated that the superfiuid-solid ^He interface 
can be set in motion with a sound wave which transmits 
through the interface, giving rise to the processes of crys- 
tallization and melting. 



SPHERICAL GEOMETRY 

In experiment a solid nucleates at the needle point 
and then grows free. A ratio of the crystal sizes in the 
different directions for the crystals grown is not drastic. 
The ratio of the maximum to minimum size does not 
exceed 2-3, see also Figs. [1] and [2] Here we consider the 
stability of the spherical shape of a growing solid. For 
simplicity, we assume the surface tension a is isotropic. 
In addition, we neglect the acceleration of gravity and, 
as above, treat the equations linearized in the interface 
perturbation. 

The equation specifying the interface motion is taken 
as r = Rsit, n) = R{t) + ({t, Q), where R{t) is the ra- 
dius of the unperturbed spherical interface and C(i, ^) = 
J2i Ci{'t)Yi{^) is an interface perturbation expanded in 
the spherical harmonics of degree / = 0, 1, 2 . . . We look 
for velocity potentials of the normal Vn = \/(f)n and su- 
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perfluid = V0s motions of the form 

Using the same boundary conditions as above, we find 
the velocities of the unperturbed flow of the liquid phase 



-R{p' ~ Ps)/Ps , 



i?, 



and coefhcients Ai (t) and Bi (t) describing the perturbed 
motion of the interface 

Employing the same boundary relation for the pressures 
in the phases and the same dependence between growth 
rate and chemical potential difference, we obtain for the 
undisturbed growth of the solid phase after some algebra 



AP 



2a 



p'-p R 



The growth equation looks exactly like the motion of 
a particle with the effective mass M{R) = AirpdR'^, 
drag force 4ttR^p'K~^R, and potential energy U{R) = 
AiraR^ - {{p' - p) / p)AP{AiiR^/?,), AP = Poc - Pc being 
overpressure. 

The small amplitude perturbations for the spherical 
surface of the solid phase are described by the relation 



Pcf 



(l^l)(l + 2) l-l , ^ 

a- '-^ p,fy-^R]Ci = 0. (11) 



i?2 



l + l 



For large radius R and Z oo so that q = l/R is fixed, 
equation pT|) goes over to the result for the plane geom- 
etry. Under steady conditions R{t) = R the spectrum of 
crystallization waves on the spherical interface is deter- 
mined from 



lUJl 



l + l 



PcfK R 



PctR^ 



(/-l)(/ + l)(/ + 2) = 0. 



Equation dTTI) can have unstable solutions giving rise to 
the generation of crystallization waves at the interface. 
The answer whether the spherical interface will be stable 
or unstable depends strongly on the history of the time 
behavior R{t). To illustrate, we consider the Richtmyer- 
Meshkov situation when a solid with radius R and inter- 
face rate R is subjected to a spherical shock acceleration 
R = V{0)S{t). In the linear regime the perturbation am- 
plitude can be estimated according to 



U(0) (l-l) 



3R + l)p'/{p,iK) 



1 



'3ft I p'(i + l) ^ 



0(0). 



Here 0(0) is the initial perturbation amplitude of l-th 
harmonic at the interface immediately after the shock 
acceleration and U(0) is an additional velocity acquired 
by the interface. The mode with Z = 1 is obviously not 
involved because this harmonic describes the displace- 
ment of a sphere as a whole. Like the case of the planar 
geometry, the initial growth of the interface perturbation 
is a linear function of time 

ci{t) = tv{m-m{^)iR- 

The long-time behavior of the interface perturbation is 
strongly governed by the magnitude and the sign of the 
growth rate R. In fact, provided that a solid '^He globe 
either grows at any rate or melts not so fast at the mo- 
ment of shock acceleration, i.e.. 



R > 



3 pciK 



the growth of the perturbation amplitude saturates. For 
larger harmonics, the stabilizing role of the finite damp- 
ing of crystallization waves increases. In the opposite 
regime when a ''He solid melts sufficiently fast, the lin- 
ear growth of the interface amplitude in time will cross 
over into an exponential increment of the perturbation 
amplitude. 

The effects in the spherical geometry can be studied 
by focusing a high-intensity sound wave in the middle of 
an experimental cell. The experiments (3l| on nucleation 
of solid *He with two hemispherical piezoelectric trans- 
ducers glued together to make a spherical geometry have 
shown that it is possible to achieve very high pressure 
amplitudes, more than 100 bar in the bulk liquid "^He. 



SUMMARY 

We have investigated here the analogue of the Ray- 
leigh- Taylor and Richtmyer-Meshkov instabilities of crys- 
tallization waves at the accelerated superfluid-solid ''He 
interface. Our analysis, made within the linear the- 
ory in perturbation, shows that the Rayleigh- Taylor and 
Richtmyer-Meshkov instabilities can occur as well as the 
parametric Faraday instability [l8| which we have treated 
as a periodically driven version of the Taylor instability. 
The plane and spherical interfaces are considered, and 
the first-order linearized equations are found for the per- 
turbation amplitudes. 

Unfortunately, it is difficult to make a well-founded 
conclusion in favor of destructing the crystal faceting as 
a result of the impulsively accelerated interface. For a 
quantitative analysis, it is necessary to have the images 
of crystal evolution taken successively in time. How- 
ever, such aspects as the Taylor-like instability cannot 
be rejected for the highly accelerated superfluid-solid ''He 
boundaries. 
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Regarding the well-faceted and atomically smooth 
crystal surfaces which may have an infinitely large surface 
stiffness, we can make the following remarks. The large 
value of surface stiffness 7 is a factor, first, for nonzero 
wave vectors q, which prevents the development of the 
Taylor instability. The most distinctive feature of the 
smooth crystal surface from the rough one is the exis- 
tence of a non-analytic cusp-like behavior in the angle 
dependence for the surface tension, e.g., 0, The 
presence of a singularity leads to qualitative distinctions 
in the development of instabilities 16|. First and fore- 



most, the threshold magnitude and the conditions deter- 
mining the development of the Taylor instability prove 
to be dependent on the initial amplitude of the interfa- 
cial perturbation. The smaller the initial perturbation 
amplitude, the larger the necessary threshold magnitude 
of acceleration, approximately as 1/Cq(0). We note here 
that the onset of the Taylor instability is favored at the 
vicinal surfaces whose orientations are tilted by a small 
angle with respect to the high-symmetry faceted ones. 
The larger the slope of the vicinal plane, the smaller the 
initial amplitude of interfacial perturbations is required 
for the development of the instability at the same mag- 
nitudes of acceleration. 

The constant acceleration of an interface works like 
an effective gravity. Hence, the physical picture for the 
Taylor case of constant acceleration is analogous to the 
gravity-driven fingering of a crystal atop a fluid and 
can qualitatively be interpreted in terms of an effective 
amplitude-dependent stiffness 7 ^ l/'zCq [16]. For suffi- 
ciently large perturbation amplitudes Cgj the difference 
between the cases of smooth and rough crystal surfaces 
disappears. 

In the case of impulusively accelerated smooth inter- 
face, the very initial stage of instability will be similar 
to that of a rough surface with the perturbation am- 
plitude increasing linearly in time. Insensitivity to the 
surface state results from the fact that the inertial prop- 
erties of the interface, associated with its effective den- 
sity, are mainly responsible for a linear response on in- 
stant loading. However, the specific time when the lin- 
ear time growth of perturbations breaks down becomes 
amplitude-dependent for the smooth surfaces and shorter 
as compared to the rough surfaces. 

Experimental and theoretical study of these effects can 
be useful for clarifying physical aspects of the crystal ^He 
growth under high drives and fast dynamics. 
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